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Based on the competition between members of a hierarchy of length scales in complex multi-scale 
systems, it is shown how clustering of active quantities into concentrated sets, like bubbles in a 
Swiss cheese, is a generic property that dominates the intermittent structure. The halo-like surfaces 
of these clusters have scaling exponents lower than that of their kernels, which can be as high as 
the domain dimension. Examples include spots in fluid turbulence and droplets in spin-glasses. 
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It has long been recognized that active quantities in 
complex systems of many types are not distributed evenly 
across a domain but cluster strongly into irregular bub- 
bles, as in a Swiss cheese. The nomenclature, the nature 
and shape of the bubbles, and the physics in each sub- 
ject is substantially different: spottiness in high Reynolds 
number fluid turbulence [1-8] and boundary layers [9]; 
droplet formation in spin-glasses [10-12]; clustering be- 
haviour in networks [13, 14]; the preferential concentra- 
tion of inertial particles [15-19] with applications to rain 
initiation by cloud turbulence [20]; and magnetic bub- 
bles in astro-physics [21], are just a few examples. These 
clustered sets display strong features whose typical length 
scales are much shorter than their averages, thus raising 
the question of the nature of the interface between them 
and their surrounding longer scale regions. For instance, 
in spin glasses the ‘surface’ of the droplets has a fractal- 
like structure whereas the droplets themselves take the 
full domain dimension [12]. In fluid turbulence the con- 
centrated sets on which vorticity accumulates are tubes 
and sheets, although the precise fractal nature of these 
is unclear. These sets dominate the associated Fourier 
spectra which display a spikiness that is the hallmark of 
what is usually referred to as intermittency [1-5]. The 
ubiquity of this irregular bubble-like topology suggests 
the existence of a set of organizing principles behind its 
occurrence in complex multi-scale systems. Using simple 
but broadly applicable mathematical ideas, this paper will 
demonstrate that the dominant physical principle behind 
clustering is the existence of a hierarchy of length scales 
whose members are in competition. 


Consider a d-dimensional system whose smallest char- 
acteristic (integral) scale L is such that the system is sta- 
tistically homogeneous on boxes Q = [0, L]. Moreover, 
it is endowed with the following two properties. Firstly, 
at each point x € Q, it possesses an ordered set of length 


FIG. 1: An illustrative slice through Q for one value of n: 
the black kernels (bad regions) are surrounded by red halos. 
Very small-scale behaviour concentrates on the black and red 
regions which constitute the set Af (Lnkn > 1). The halos 
have scaling exponents lower than those of the black kernels. 


scales £,, = ln(x) associated with a hierarchy of features 
labelled by n > 2 
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The £, could be thought of as an ordered set of correla- 
tion or coherence lengths; their inverses k,,(x) = & (x) 
clearly obey 1 < Lkn < Lkn41. The second assumption 
is that the ensemble averages of the Lkn(x) are bounded 
above by some ordered, positive parameters of the system 
satisfying 1 < Rn < Rn+1 


1< L(kn) < Ran. (2) 


The ensemble average (-) is a spatial average with re- 
spect to the Lebesgue measure over Q. Thus, while the 
ordering of the £ (x) must be respected at each point, 
the Zn themselves could be quite rough; e.g. they could 
consist of a series of step functions. If they become very 
small near points x* then they must obey £», > O(r¢~*) 
(r = |x — 2*| and £ > 0) so as not to violate (2). 


Following an idea used in ref. [22], consider the real 
arbitrary parameters 0 < u < 1 and 0 < a < 1 such that 
puta = 1. Use Hölder’s inequality, |AB| < ae hee 
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where 5+ 7 = 1, with p = 47 and q=a_ 
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Re-arranging and factoring out a term (K°) gives 
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Lower bounds on the ratio (K2) / (Kn) can be found from 
(2) thereby turning (4) into 


(=y E [eten >0. (5) 


While it is possible that the integrand in (5) could 
be positive everywhere in Q, this cannot be assumed; 
the generic case is that the integrand could take either 
sign[35]. With the definition £L, = L Rp “” we have the 
pair of inequalities 
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for which > is valid on regions where the integrand is 
positive, designated as good regions, and negative (<) on 
bad regions. The term (LnkKn)” on the right hand side of 
(6) remarkably contains the arbitrary parameter u which 
lies in the range 0 < u < 1. Its existence is important 
because the ordering in (1) makes it clear from (6) that 
everywhere within the bad regions (<) there are large 
lower bounds on kn with exponents containing 1/1 
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Let A} be the set on which Lnn > 1 and A; the set 
on which Lakn < 1. Then all the bad regions (<), des- 
ignated by the clusters of black kernels in Figure 1, lie 
in A}. The red halos also lie in A}, and correspond to 
those parts of the good regions (>) neighbouring the bad. 
It is in these halos where the lower bound (LnkKn)” be- 
comes operative. The white areas of Figure 1 correspond 
to A, in which the «kn can be randomly distributed sub- 
ject to their ordering in (1). It is clear from (5) that the 
existence and location of the clusters may differ for each 
n. A physical picture that displays all clusters for every 
n would be the union At = UA. 

To show that the volume V+ of A} comprises a small 
part of Q, Chebychev’s inequality relates the normalized 
Lebesgue measure m( A7) to the integral of Lk, over Aj 


f Lkn dm > m( A$) RY" = LVE RY”. (8) 
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Together with the relation ve Lrndm < (Lkn) < Rn 
we have 


m(At) < RP". (9) 


Hence m( A} ) is significantly smaller than unity and de- 
creases as R, increases. Thus A} can fill, at most, a 
small fraction of 2. With such sparse information it is 
difficult to estimate the Hausdorff or the fractal dimen- 
sions of A}, but it is still possible to estimate scaling 
exponents [23]. This entails making a third assumption 
of self-similarity to estimate the smallest number of balls 
N+ of radius AF needed to cover A}. Defining AF as 


(At)? = ke = (Why /?, (10) 
for some p > 1, it is clear that kf cannot be large enough 
when p = 1 because of (2). However, any value[36] of 
p > 1 will do that makes kt large enough to be a member 
of A}. The simplest and worst estimate would be to 
write 


Naw (E/N). (11) 


Inequality (9), however, shows that Af occupies only a 
small fraction of Q. A multiplicative factor of m( A} ) is 
introduced thus 
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i; = m(At)(LkE)*. (12) 


Instead of using (9) to estimate m( A} ), an assumption 
of self-similar scaling is introduced that requires that the 
change in volume of the balls with respect to n should 
scale as V+ (the volume of Af) scales to L¢. Thus 


d 
m(At) A (= , (13) 


We observe that the definition of the set A; in principle 
involves the length scales L and A}, but not overtly Af}1- 
Yet the good and bad sets involve all three scales; L, 
AF and A*,,. The self-similarity assumption (13) is an 
assumption about the nature of the set A;* that relates 
successive length scales A and Af}; in an ad hoc, yet 
reasonable, fashion. Using (13) in (12) we have 
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From these, two estimates for Mọ emerge, one each for 
the red halo and black kernel regions of Figure 1, whose 
scaling exponents[37] are independent of p 
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TABLE I: Summary of conclusions regarding the sets A; and 
the coloured regions in Figure 1. 


Figure 1 black red white 
Set An An An 
Inequality (6) < (bad) > (good) > (good) 
Exponent <d <d(1— u) =d 


For the former, the > direction of the inequality in (6) 
has been used together with a simple Hölder inequality 
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whereas for the latter kpn < Kn+ı has been used. In 
contrast, without any evidence of contraction of volume, 
the formula corresponding to (12) for Mz is 


New ei = (Lrk,)° RY, (17) 


where k, satisfies £,4, <1. The uniform scaling expo- 
nents in (15) are bounded by 
Dt 
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n,ker < d (18) 
whereas D` = d from (17). The coefficients R” in (14) 
to (17) reflect the fact that this effect is taking place 


only at length scales smaller than LR, bi t. The red halo 
clearly plays the role of an interface of small but finite 
thickness between the d-dimensional (white) outer region 
and the (black) inner kernel whose dimension can be as 
high as d but could be less. When Dr ker Saturates its 
upper bound we have 
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< d(1 => u) < Dig =d. (19) 
For the red region to have an exponent at least d — 1 
(a surface), u would lie in the range 0 < u < 1/d. 
Without equations of motion, a numerical experiment 
would be necessary to estimate the Rn by finding the 
maximum value of the ensemble average (Kn). In prin- 
ciple u could then be found from numerical estimates of 


gerit ~ LRn 1/4 within the black kernels although if the 
Kn take very large values there it might not be possible 
to achieve resolution. p itself may have upper and lower 
bounds that are themselves n-dependent, as in ref. [22]. 

We now proceed to discuss some examples. The first 
ideas on clustering came more than half a century ago 
from Batchelor and Townsend [2] who observed intermit- 
tent behaviour in their high Reynolds number flow exper- 
iments, closely followed by observations in boundary lay- 
ers by Emmons [9]. Batchelor and Townsend called this 
phenomenon ‘spottiness’ and suggested that the energy 
associated with the small scale components is distributed 
unevenly in space and roughly confined to regions which 
concomitantly become smaller with eddy size [24]. Man- 
delbrot then suggested that these clustered sets on which 


energy dissipation is the greatest might be fractal in na- 
ture [25]. In measurements of the energy dissipation rate 
in the atmospheric surface layer, Meneveau and Sreeni- 
vasan interpreted the intermittent nature of their signals 
in terms of multi-fractals [3]; a newer generation of exper- 
iments measuring intense dissipation in turbulent flows 
have been pursued by Zeff et al [4]. Sreenivasan and 
Bershadskii have recently suggested that the clustering 
of high frequencies in a turbulent signal can be charac- 
terized by a scaling exponent [5]. 

The extremely rapid time evolution of sets of high vor- 
ticity or strain in fluid turbulence is an important is- 
sue; many computations exist showing how these take 
on the nature of quasi-one-dimensional tubes and quasi- 
two-dimensional sheets which have short lifetimes [6, 8]. 
An alternative to studying the problem in a statistical 
manner is to include time in the ensemble average (-), in 
which case the semi-infinite nature of the time-axis sug- 
gests a different measure might be necessary[38]. With 
specific reference to the Navier-Stokes equations, analysis 
is not advanced enough to deal with the full space-time 
equations (except see ref. [26]); conventional methods of 
analysis use Sobolev norms to L?-average over space and 
remove the pressure [6, 27, 28] leaving only time. In ref. 
[22] a hierarchy of kn have been constructed which are 
comprised of ratios of norms (of derivatives of order n) 
and therefore functions of time only; thus the clusters 
of Figure 1 are merely gaps in the time-axis. It is then 
necessary to prove that they are finite in width and de- 
creasing with increasing Reynolds number. This involves 
finding bounds on p. 

The second example is that of the low-temperature 
phase of spin glasses [29, 30]. Our conclusions regard- 
ing the fundamental role played by the competition be- 
tween members of a hierarchy of length scales is con- 
sistent with the observation of ultrametricity in spin- 
glasses, a term that is used to denote the presence of 
a hierarchy of scales [30, 31]. This has been observed in 
computations on the low-temperature spin glass phases 
of the Sherrington-Kirkpatrick [29, 31, 32] and Edwards- 
Anderson models [33], as well as in dynamic phenomena 
in complexity [34]. The results of this paper, and par- 
ticularly with reference to (19), are consistent with the 
droplet theory [10-12] where the kernel of the droplet is 
of full dimension d but its surrounding ‘surface’ has a 
scaling exponent < d. In fact, Palassini and Young [12] 
have shown numerically that Diis = 2.58 + 0.02 when 
d = 3 and DË jo = 2-77 0.02 when d = 4. 


In conclusion, we have shown that in a system endowed 
with a competitive hierarchy of correlation lengths, a 
clustering effect ensues in which length scales smaller 
than a critical value, and much smaller that the ensem- 
ble average scale, aggregate into small intense regions. 
The kernels of these intense regions are surrounded by 
halos that have scaling exponents smaller that of the do- 
main dimension d. We have expectations that this idea 
of competition between scales may be a useful paradigm 
in explaining the behaviour of multi-scale systems. 
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